C ERN-TH/98-69 
gr-qc/9804076 



Seeds of large-scale anisotropy in string cosmology 

R. DurrerW, M. Gasperini^ 2 ) , M. Sakellariadou^ and G. Veneziano^ 3 ) 
Departement de Physique Theorique, Universite de Geneve, 

24 quai E. Ansermet, CH-1211 Geneva, Switzerland 
^ Dipartimento di Fisica Teorica, Universita di Torino, 
Via P. Giuria 1, 10125 Turin, Italy 
(3) Theory Division, CERN, CH-1211 Geneva 23, Switzerland 
' (April 1998) 



5h 
Oh! 



Pre-big bang cosmology predicts tiny first-order dilaton and metric perturbations at very large 
scales. Here we discuss the possibility that other - more copiously generated - perturbations may 
act, at second order, as scalar seeds of large-scale structure and CMB anisotropics. We study, in 
particular, the cases of electromagnetic and axionic seeds. We compute the stochastic fluctuations of 
their energy-momentum tensor and determine the resulting contributions to the multipole expansion 
QQ ' °f the temperature anisotropy. In the axion case it is possible to obtain a flat or slightly tilted blue 

£SJ ' spectrum that fits present data consistently, both for massless and for massive (but very light) 

axions. 



I. INTRODUCTION 



> 
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String theory has recently motivated the study of a cosmological scenario in which the universe, starting from the 
string perturbative vacuum, evolves through an early inflationary "pre-big bang" phase |Q, until a transition to the 
1 radiation-dominated, decelerated evolution occurs. 

In spite of some attractive aspects of the pre-big bang picture, such as the underlying duality symmetry which 
naturally selects perturbative initial conditions and automatically leads to inflation [|l]||], it is fair to say that such 
O" 1 a cosmological scenario is far from being understood in all of its aspects. In particular, on the more theoretical 
, side, one is lacking a complete and consistent description of the high-curvature, strong coupling regime, where the 
transition from the pre- to the post-big bang era is expected to take place [Q. Furthermore, opinions vary as to 
whether or not the pre-big bang scenario needs a large amount of fine-tuning. On a more phenomenological side, the 
main outstanding problem is to reproduce the observed amplitude and slope of the large-scale temperature anisotropy 
H and of large-scale density perturbations. The difficulty is that, unlike in the more conventional (de-Sitter-like) 
inflationary picture, the amplification of scalar and tensor metric perturbations here leads to primordial spectra that 
grow with frequency Q, and whose energy density is normalized to an almost critical value at some short scale Q 
(typically the GHz); in this way, too little power is left at scales that are relevant for anisotropics in the Cosmic 
Microwave Background (CMB) || or to the problem of large-scale structure (unless the high-curvature phase is long 
enough and characterized by an almost constant dilaton field j)| ) . 

In this paper we address this problem and we discuss a possible solution, based on the contribution of "seeds" |lC| ] 
to density fluctuations and to the anisotropy of the CMB radiation. The seeds are produced, in our context, by the 
amplification of quantum fluctuations of some other fields, which are present in string theory, but are not part of the 
homogeneous background whose perturbations we wish to study. 

We shall consider two examples, in which the seed inhomogeneity spectrum is due, respectively, to vacuum fluctu- 
ations of the electromagnetic (EM) |ll|] and of the (Kalb-Ramond) axion (AX) |l2| field. Both cases are typical of 
string cosmology, since no inhomogeneity is produced, in either case, in a conventional scenario based on Einstein's 
equations, without axion and dilaton. The spectra of EM and AX perturbations can be much flatter than those of 
scalar and tensor perturbations of the metric and of the dilaton field. 

The idea of using the EM fluctuations as seeds was already discussed in a previous paper jljj, using however the 
perfect fluid approximation for the EM stress tensor. Here we will compute the scalar components of the energy- 
momentum-tensor fluctuations due to the EM and AX seeds including an important anisotropic stress term, and 
will relate them to the primordial spectral energy distributions. When these seed inhomogeneities are inserted in 
the perturbed Einstein equations they generate scalar-metric fluctuations which are largely controlled, for seeds with 



1 



small enough anisotropic stresses, on super-horizon scales, by the so-called compensation mechanism [H. Finally, 
scalar-metric perturbations can be converted in a standard manner into temperature fluctuations AT/T via the Sachs- 



Wolfe effect 15 1 . We will discuss whether the metric perturbation spectrum induced by seeds can be flat enough to 
match present observations, consistently with the COBE normalization of the amplitude on large scales, and with the 
high-frequency normalization of the primordial seed spectrum. 

The paper aims at being rather self-contained and readable by non-specialists in string and/or cosmological pertur- 
bation theory, and is organized as follows. In Section || we set up the relevant equations needed to study super-horizon 
perturbations in the presence of seeds, and give their generic solution for seeds with "small" or "large" anisotropic 
stresses. We also discuss the way the perturbations enter the multipole expansion of AT/T via the Sachs- Wolfe effect. 



In Section III, after recalling known results about scalar, tensor, electromagnetic and axion perturbations in the pre- 



big bang scenario, we estimate the contribution of the two latter sources to the fluctuations of the energy-momentum 



tensor, including the case of massive axions. In Section IV we combine the results of the previous two sections and 



compute the contribution of EM and AX seeds to AT/T. Using COBE data, we finally discuss, in the various cases, 
whether the seed mechanism alone is able to give a satisfactory explanation of large-scale temperature anisotropies. 
Section contains our conclusions. Some technical details are relegated to the three appendices. 

Notation: The Friedmann metric is given by a 2 (—dri 2 + j i jdx l dx : '), where a denotes the scale factor and r\ is con- 
formal time. Spatial indices, 1, 2, 3 are denoted by latin letters while spacetime indices, 0, 1, 2, 3 are denoted by greek 
letters. A dot denotes derivative with respect to r\. 

II. LARGE-SCALE PERTURBATIONS IN THE PRESENCE OF SEEDS 

Before calculating CMB anisotropies for specific examples in the context of string cosmology, we derive a general 
formula for large-scale CMB anisotropies in models with seed perturbations. 

A. Cosmological Perturbation Theory with Seeds 

In this subsection we give a brief reminder of gauge-invariant perturbation theory with seeds. More details can be 
found in Refs. JTcjJT^ ] . By seeds we mean an inhomogeneously distributed form of energy, which contributes only a 
small fraction to the total energy density of the universe and can thus be considered as a perturbation. Furthermore, 
we consider seeds that interact only gravitationally with the cosmic fluid. 

We restrict our discussion to scalar perturbations, which arc of primary interest here. The corresponding equations 
for vector and tensor perturbations can be found in |l6|] . The metric of a perturbed Friedmann universe is 

where g(°) denotes the unperturbed metric: 

g^dx^dx" = a 2 (r])(-drf + ltj dx l dx j ) . (2.2) 

Here a is the scale factor, r\ denotes conformal time and 7 represents a metric of constant curvature K — ±1, 0. Since 
we will be interested in a Friedmann universe that has undergone substantial inflation, we neglect K in the sequel, 
setting 7y = 

For scalar perturbations, a Fourier component of h^ v with wave vector k can by parametrized by 4 scalar functions 
A, B, Hl and Ht, defined by 

h(k) = h^{k)dx^dx u = -2A(fc)(d7y) 2 - 2i^-B{k)dj 1 dx 1 



H L (k) + l -H T {k) 



dud^dx 1 - 2^fy-H T (k)dx l dx j . (2.3) 
k 



These four functions are gauge-dependent, i.e. they depend on the choice of coordinates. In order to define gauge- 
independent metric variables, we first make use of two geometric quantities: the spatial part of the scalar curvature of 
the perturbed metric, SR, and the shear (traceless) part of the extrinsic curvature, x^ anlso \ An elementary calculation 
gives fu|: 

5R = Ak 2 a - 2 n , K = H L + \h t , (2.4) 
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4"" SO) = ak - l^J a, a = H T /k-B. (2.5) 

Studying the gauge transformation properties of A, 1Z, and a, one easily finds that the following variables, called the 
(Fourier components of the) Bardeen potentials, are gauge-invariant (see [ ]l7| , |l8| | ) : 

$ = 11 - {a/a)k- l a , (2.6) 



* = A- (a/ajk^a-k^a 



(2.7) 



(Note that, throughout this paper, we shall always express the Bardeen potentials in momentum space, even without 
indicating their k dependence explicitly.) 

Next, we discuss the perturbations of the energy-momentum tensor. Let us define the perturbed energy density 
pKpert) an( j 4_ ve i oc ity field u as the time-like eigenvalue and eigenvector of the energy-momentum tensor: 



u z = -1 



The Fourier components of the perturbations in the density and velocity field are determined by 

p ( P ert) = p{1 + J) ^ 



(2.8) 

(2.9) 



u° = (1 - A) , 



it 



(2.10) 



where p denotes the unperturbed background density. The temporal component vP is fixed by the normalization 
condition. We project the stress tensor onto the 3-space orthogonal to u: 



up p °t 



P„ 



and define the scalar perturbations of r by: 



P 



1 +TTL + TtTTT I S. 



k,w 



(2.11) 



(2.12) 



The variable ttl describes the pressure perturbation, ttt is the potential of the anisotropic stresses and p is the unper- 
turbed background pressure. Studying the behaviour of the quantities 5, v, ttl and ttt under gauge transformations 
|n_9[, one finds the gauge- invariant variables: 



n 



7TT, 



r = ttl - {c 2 Jw)5, V = v-k- 1 H 7 



D = 5 + 3(1 + w)(a/a)k- x (y + a), D g = 5 + 3(1 + w)TZ 



(2.13) 



Here n is the anisotropic stress potential, T is the entropy perturbation, V is the peculiar velocity potential, D and 
D g arc different choices for a gauge-invariant density perturbation variable (for a physical interpretation of these 
variables, see 18 ; l|]). Finally, w = p/p denotes the enthalpy and c 2 s = p/p stands for the adiabatic speed of sound. 
In this paper we shall limit ourselves to adiabatic perturbations (r = 0). 

The perturbation of Einstein's equations and of energ y-mom entum conservation can be expressed in terms of these 

. We obtain two constraint equations: 



gauge-invariant variables (a derivation can be found in 18 

AuGa 2 pD = fc 2 $, 
AuGa 2 {p + p)V = k [(a/a)* - $ 



(2.14) 
(2.15) 



two dynamical equations: 



-8nGa 2 pIl = k 2 ($ + *), 

&7TGa 2 p [r + {c 2 Jw)D g + (2/3)fc 2 n] = 



(2.16) 



and two conservation equations: 



D a — 3w a (a/a)D a — —k [(1 + w a )V a + 2(d/ 'a)w a k 1 H a \ 
+3(l + w a )4irGa 2 (p + p)(V -V a ) , 



V a + (a/a)V a 



1 + w a 



kD a 



1 + w a 



— kT a + k^ - 



2w n 



3(1 + w a ) 



kU a 



(2.18) 
(2.19) 



The above conservation equations hold for any component a of the fluid stress-energy tensor which interacts with 
the other components of the cosmic fluid only gravitationally. The variables c a and w a denote the adiabatic speed of 
sound and the enthalpy of the fluid component, respectively. The total perturbations are defined as the sums: 



pD = ^p a D a , (p +p)V = ^(/? a +p a )V a , etc. 



(2.20) 



For interacting matter, the corresponding equations can be found in |lq ]. 

In order to complete the above analysis we also need equations of state for the matter sources, which relate for 
instance T and II to D and V. Due t o the Bianchi identities, the conservation equations for the total cosmic fluid 
follow from the fi eld e quations ( 2.14 )- (|2.17 ). Thus, we need n ot m ak e exp licit use of both dynamical equations, but 
we can use, say, (2.16) and one of the conservation equations ( [2.18 ), ( 2.1 9| ) for the total fluid. 

We now add to the perturbation equations an inhomogeneous energy-momentum distribution, Tfi), generated by 
seed fields that do not interact with the cosmic fluid other than gravitationally. 

Since, by definition, seeds do not contribute as sources of the homogeneous background, the energy-momentum 
tensor T p s J is gauge- invariant by itself pc[| , and can be calculated by solving the field equations for the seeds in the 

( s) 

unperturbed background geometry. Let us assume that we can express the Fourier components of T pv in terms of four 
scalar "seed-functions" f p , f p , f v and / ff (we just neglect vector and tensor contributions; since they are decoupled 
from density perturbations, in the linear approximation, this will not affect our results for scalar perturbations): 



li s) (k, )) )=ay = M%(k ))) ), 



— n 2 n( s ) - 



L j0 1 



M) 



-iM%f v (\<L,r)) , 



= M J 



A*fc* n ( s ) 



f p (k, t}) + y / ff (k, 77) ) jij - kikjUik, 77) 



(2.21) 
(2.22) 



(2.23) 



Note that f p and f p have dimension £~ 2 , while /„ has dimension £ _1 and is dimensionless. Here M denotes an 
arbitrary mass scale, introduced for dimensional reasons, which will eventually drop out in physical predictions. 

Given an energy- momentum tensor T pv , which in general contains vector and tensor contributions, the scalar parts 
/„ and /,r are determined by the identities: 



ik j T$ = M 2 k 2 f v , 



l 0j 

2„,klrp(s) 
kl 



--M 2 k 4 f., 



(2.24) 



On the other hand, f v and / ff are related to f p and f p , by the conservation equations V v T p s J = 



f P + k 2 f v + (h/a)(f p + 3f p ) = Q, 



(2.25) 



/„ + 2(d/o)/„ - f P + (2/3)fc 2 / 7 r = 



(2.26) 



In the presence of seeds, and in the approximation in which perturbations are treated linearly, the total geometric 
perturbations can be separated into a part induced by the seed s, $ s , a nd a part induced by the perturbations of 



the cosmic fluid, 9 m ,$> r 



The perturbed Einstein's equations ( 2.14 ) and ( 2.16 ) become 

fc 2 $ = 4TrGpa 2 D + e[f p + 3(d/a)/„] , 
$ + * = -8irGa 2 k- 2 P Il - 2e/ w , 



(2.27) 
(2.28) 
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where e = 4irGM 2 . If we define 



9 r . 



(2.29) 



with: 



we easily find 



k 2 d> s =e[f p + 3(d/o)/„)] 



= 4nGpak 



D g + 3{l + w) 



V 



3(l+io)$ 



-$ m - 87rGa 2 pIir 



(2.30) 

(2.31) 
(2.32) 



Equation (2.31) has been written in terms of the gauge- invariant density perturbation D g , because this choice will 
simplify our final equations. Physically, D g corresponds to the density perturbation in the flat slicing. The evolution 
of D g and V is described by the conservation equations ( 2.1 8| ) and ( 2.19| ), which read explicitly: 

-(l + w)kV, 



D g + 3(c 2 



W)-Dg 

a 



V A- -(1 - 3c 2 s )V = jfe(* - 3c 2 $) 



1 



-D„ 



2w 



3(1 



-MI 



(2.33) 
(2.34) 



To simplify the analysis, we will assume w = c 2 = constant. The unperturbed background equations are then solved 
by a oc ?y r , with r = 2/(3w + 1). Since we are interested in very large scale perturbations in the c osmic micr owave 
background, we concentrate our discussion on super-horizon scales, such that kr\ <C 1. Eqs. (2.31) and (2.29) then 
lead to 



$ = 



1 



3(1 



w 



9r 2 (l 



w 



-{kri) 2 ^ s 



(2.35) 



where r — 1 for the radiation-dominated era, and r = 2 for the matter-dominated era. 

The evolution equation for D g , Eq. ( [2.33 ), implies dD g /d(kri) = —(1 + w)V. In the physical picture we have 
in mind, metric perturbations are triggered by the presence of the seeds alone, and we do not want to include an 
arbitrary contribution from the perturbations of the ho mogen eous sources. We thus require D g (0) — 0, which implies 
D g ~ kr/V. Hence, we m ay negl ec t the D ff -term in Eq. ( 2.35| ) for krj <gC 1. 

Combining Eqs. (fE3§), ( gjg ), ( gjgg ) we find, on super-horizon scales, 



dV r 2w s , 2(fcrj) 2 

\T/ = 1 VA (kn) $ A - — 

d(kr)) k V 3r 2 (l + w) V /; 9r 2 (l + w) 



(2e/ w + * + $) 



(2.36) 



The two equations ( |2.35| ) and (2.36) relate the three variables ty, $ and V once the seeds are given. To proceed, 
we need an equation of state to close the system. For single component fluids this equation usually takes the form 
II = H(D g ,V). We are interested in large-scale CMB anisotropies, which are induced at recombination and later, 
when the universe is already matter-dominated, with p <C p- Thus, in what follows, we will consider the case II = 0, 
which implies 



$ + * 



(2.37) 



Furthermore, in a matter-dominated Friedmann universe, r — 2 and w = 0. The equation of motion for V, obtained 
by combining Eqs. (p.35|), (2.36), (2.37), then reads 



^ + ±V = -±(k V )H s - 2eU = [fp A- 3(a/a)f v] 



2e/ w 



(2.38) 



In the next subsec tion^ w e sha ll see that t he large-scale anisotropies of the CMB are determined by the combination 
Using Eqs. fl2.35[), (pT37j) and (^.38[), we find immediately: 



-^-^ 2 e[f P + 3(a/a)f v ] 



(2.39) 
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Modulo numbers of order unity, which can be computed case by case, we finally arrive at the estimates: 



dV V f , / a 



(2.40) 



Depending on whether ij 2 (f p + 3(a/a)f v ) or dominates in Eq. ( 2.40 ), we can distinguish between seeds with small 
and large anisotropic stresses. We will discuss in Section III to which case o ur str in g-cos mology seeds belong. 
If the term erf (f p + 3(a/a)f v ) = x 2 $ s dominates, we conclude from Eqs. ( 2.38| ),( 2.39 ) that 



(2.41) 



on super- horizon scales. This suppression of the total geometric perturbations, if compared with th e sou rce pertur- 
bations alone, is known under the name of "compensation" iQ. The conservation equations (2.18), (2. IE) show that 
the presence of seeds induces matter perturbations that try to compensate the gravitational potential of the seeds. 
Since anisotropic stresses in the seeds cannot be compensated by a perfect fluid, compensation is not effective, if 
anisotropic stresses dominate. But, as shown here (see also [fblf), the phenomenon of compensation is quite generic 
and, to a large extent, independent of the spectrum of seed fluctuations. 



B. The Seed Contribution to CMB anisotropics 



In this subsection we calculate the CMB anisotropics for models where perturbations are induced by seeds, and 
their contribution to AT/T via the Sachs- Wolfe effect |15 . We first discuss in general the motion of photons in a 
perturbed Friedmann universe. 

We make use of the fact that the equations of motion of photons are conformally invariant. More precisely, two 
metrics that are conformally equivalent, 



ds 2 = a 2 ds 2 



(2.42) 



have the same light-like geodesies, only the corresponding affine parameters are different. Let us denote the two affinc 
parameters by A and A respectively, and the tangent vectors to the geodesic by 



dx 



dx 



n= d\' n = Tx 



n 2 = n 2 = , n° 



1, n 2 



1. 



Setting n = 1 + 5n° , the geodesic equation for the perturbed metric 



ds 2 



{rjftu + h tiv )dx ll dx v 



yields, to first order, 



01/ 



i/ 1 



/ 



hit^n^n^dX 



(2.43) 



(2.44) 



(2.45) 



On the other hand, the ratio of the energy of a photon measured by some observer at tf to the energy emitted at tj is 

E f _ (g ■ ")/ _ T f ( n ■ u )f 



E; 



(n ■ u)i Ti (n ■ u)i 



(2.46) 



where Uf and ii$ are the four- velocities of the observer and emitter respectively, and the factor Tf/Ti is the usual 
(unperturbed) redshift, which relates n and n. The velocity field of observer and emitter is given by 



u = (1 - A)d v + v l d t 



(2.47) 



An observer measuring a temperature To receives photons that were emitted at the time r\dec of decoupling of matter 
and radiation, at the fixed temperature Td ec . In first-order perturbation theory, we find the following relation between 
the unperturbed temperatures Tf, Ti, the measurable temperatures T Q , T^ec, an d the photon density perturbation: 



II _ _ 

Ti Td, 



Tf Ti 



Tp 

Tdec 



(2, 



G 



where 8^ is the intrinsic den sity p erturbatio n in t he radiation and we u sed oc T 4 in the last equality. Inserting 
the above equation and Eq. ( 2.45 ) into Eq. ( 2.46 ), and using Eq. (2.2) for the definition of h^ v , one finds, after 
integration by parts |n6[ : 



T 

Tdec 



\ D h) + v i™) n i + w _ $ 

4 y ■> 



(* - $)dA 



(2.49) 



Here Dg 1 ^ denotes the density perturbation in the radiation fluid, and is the peculiar velocity of the baryonic 



matter component (the emitter and observer of radiation). The final time values in the square bracket of Eq. (2.4E) 
give rise only to monopole contributions and to the dipole due to our motion with respect to the CMB, and will be 
neglected in what follow s. 

Evaluating Eq. ( 2.49| ) at final time ijq (today) and initial time r/dec, we obtain the temperature difference of photons 
coming from different directions n and n' 



AT _ ST{n) ST(n') 

rj~i rjn rji 



with temperature perturbation 



ST(n) 
T 



4 9 



Vj 'n J 



f - $ 



'(11 



(Vdec, x dec ) + / (* - x(rf))dri , 



(2.50) 



(2.51) 



Ida 



where x(ry) = xo — (rjp — r?)n is the unperturbed photon position at time rj for an observer at xq, and x^ec = ^-(Vdec)- 



The first term in Eq. (2.51) describes the intrinsic inhomogeneities on the surface of the last scattering, due to acoustic 
oscillations prior to decoupling. In general, it also contains contributions to the geometrical perturbations. This is 
especially important in the case of adiabatic inflationary models plfl . However, for perturbations induced by seeds, 
which satisfy the initial condition D g (k, 77) — > for 77 — > 0, the geometrical contributions to D g can be neglected. The 
second term describes the relative motions of emitter and observer. This is the Doppler contribution to the CMB 
anisotropics. It appears on the same angular scales as the acoustic term, and we thus call the sum of the acoustic 
and Doppler contributions "acoustic peaks" . 

The last two terms are due to the inhomogeneities in the spacetime geometry; the first contribution determines 
the change in the photon energy due to the difference of the gravitational potential at the position of emitter and 

(r) 

observer. Together with the part contained in D g they represent the "ordinary" Sachs- Wolfe effect. The second 
term accounts for red-shift or blue-shift caused by the time dependence of the gravitational field along the path of 
the photon, and represents the so-called Integrated Sachs- Wolfe (ISW) effect. The sum of the two terms is the full 
Sachs- Wolfe contribution (SW). 

On angular scales 0.1° ~ 9 ~ 2°, the main contribution to the CMB anisotropics comes from the acoustic peaks, 
while the SW effect is dominant on large angular scales. On scales smaller than about 0.1°, the anisotropies are 
damped by the finite thickness of the recombination shell, as well as by photon diffusion during recombination (Silk 
damping). Baryons and photons are very tightly coupled before recombination, and oscillate as a one-component 
fluid. During the process of decoupling, photons slowly diffuse out of over-dense regions into under-dense ones. To 
fully account for this process, one has to solve the Boltzmann equation for the photons (see, e.g. p6|). 

The angular power spectrum of CMB anisotropies is expressed in terms of the dimcnsionlcss coefficients CV, which 
appear in the expansion of the angular correlation function in terms of the Legendre polynomials Pf. 



? (n) ? (n ' , / l „„,_ * 



-£(2* + l)QP/(cos0) 



(2.52) 



Here the brackets denote spatial average, or expe ctatio n values if perturbations are quantized. 
To determine the Cg we Fourier-transform Eq. (2.51), defining 



1 



and using the identity 



e izcosl3 = ^2(2e + l)i i ji{z)P e (cose) 



(2.53) 



(2.54) 
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(where ji is the spherical Bessel function of order I). For the coefficients Ce of Eq. ( 2.52| ) we obtain: 



where 



2£+ 1 



Ct 



2 4 r) ( k >%e 



(|A.(k)| 2 
(2^+1)' 



-k z dk 



(2.55) 



(*-$)(k, %ec ) 



■ j' e (kr]o)V r (\s., r] dec ) 



+ / (* - $)(k, tj')^ - W) <V 



= |-Dg r) ( k ) ? ?<ie C )^( fc? 7o) - j£( fc7 ?o)V r (k,77 dec ) 



I/O 



(2.56) 



and j'^ stands for the derivative of j£ with respect to its argument. On large angular scales, krjdec 1 (which 
corresponds to ^ <C 100), the SW contribution dominates: 



Cf w = - I k l ,W~ 



Let us approximate the Bardeen potentials on super-horizon scales by a power-law spectrum: 

(|*-$| 2 ) = C 2 (fc) {kr,) 21 . 



(2.57) 



(2.58) 



Furthermore, we consider models where the seed contribution does not grow in time on sub-horizon scales. In this 
case the Bardeen potentials, inside the horizon, are dominated by the cold dark matter contribution, which leads to 
time- independent $ and , 3>. We can thus approximate the Bardeen potentials by 



* - $ 



C(fc)(M 7 , krj<£ 1 
C(k) , kr] > 1 



(2.59) 



We further assume that also C{k) is given by a simple power law. Thus, for dimensional reasons, it has the form 



C(k) 



_ J Mk~^ 2 (k/k 1 ) a , k< ki 



1°. 



k > ki 



(2.60) 



where TV is a dimensionless constant, and k\ denotes a comoving cutoff scale, i.e. the maximal amplified frequency 
dete rmine d by the explicit mechanism of seed production (in the case a = no cutoff is needed) . Inserting this in 
Eq. (EE57D, 



SW 



a 



where, setting x = kr], x = krj , x dec = kr, de 



N 2 - I ? I ^ 



2n 



I(k) 



dxx 1 j'i(xo — x) 



W)\ s 



dxj' e (x - x) 



dxx~<j' e (x - x) +ji(x - 1) . 



(2.61) 

(2.62) 
(2.63) 



We can see explicitely from this equation that the relevant contribution of each mode to the CMB anisotropy comes 
while the mode is still outside the horizon (krj < 1). We now distinguish two cases. 

If 7 > —1 the lower bound in eq. ( [2.63| ) can be safely extended to 0, and the integral is dominated by the region 
krj ~ 1, so that: 



H k ) ~ ji( X - Xdec) ~ ji(x ), x dec < 1 < X 



(2.64) 
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Inserting this in Eq. (2.61), the integral can be performed exactly (assuming 770 k\ £), with the result, for a < 1, 



r(2-2a) 



T(£ + a) 



4( 1 -«)r(3/2 - a) T(i + 2 - a) ' 



a < 1 



(2.65) 



(if a > 1, the integral grows towards large k and is dominated by the contributions at k ~ k\, leading to an 
^-independent result of order (N j 'kirjo) 2 ) . Comparing the above equation with the standard inflationary result p2|, 



r sw r(l-l/2 + n/2) 
1 r(^ + 5/2-n/2) ' 



(2.66) 



where n denotes the usual spectral index, we find that a is related to n by a = (n — l)/2. The scale-invariant 
spectrum, as it has been observed by the DMR experiment aboard the COBE satellite requires 



0.8 < n < 1.4 

so that, allowing for generous error bars, the COBE observations imply 

- 0.1 < a < 0.2 , 7 > -1 . 



(2.67) 



(2.68) 



Consider now the second case, 7 + 1 < 0. The integral ( 2.62 ) is now dominated by its value at the lower boundary 
and we get 



1 



( 1 + 1 y x dec |_2£+1 
If also a + 7 < 0, the fc-integral converges and we obtain (see Appendix A): 

M 2 r(-2(a + 7 )) (mec\ 2{l+1) n ,-2 a r(e + i + a + 1 ) 



2(7+1) 



i^-i(^o) 



£ + 1 

2£+l 



(2.69) 



22(«+7)( 7 + i)2r(i/ 2 _ (a + 7)) 



1 



r(£ + l-Q!-7) (2^+l) 2 



£ 2 (£-a-j) 2£(£+l)(l/2 + a + j) (I + l) 2 {l + 1 + a + 7) 



l + a + ^i 



(1/2 -a -7) 



(t+l-a-i) 



(2.70) 



Comparing again this last r esult with that of standard inflation, Eq. ( 2.66 ), and neglecting the weak ^-dependence of 
(21 + 1)" 2 [- ■ •] in Eq. ( pJ0| ), we obtain 



n - 3 + 2(q + 7) 



a + 7 < 



(2.71) 



(If, on the contrary, a + 7 > 0, the coefficients Ci are dominated by the large k (i.e. small-scale) contribution, even for 
the very low values of i. In this case the small-scale perturbations become too large, which is excluded observationally 
by the fact that the spectrum, for CMB and matter perturbations, must be close to the Harrison-Zel'dovich spectrum 

&• 

The observational limits on n thus impose 



0.1 < 7 + 1 + a < 0.2, 7 < -1 , n~3 + 2(a + 7) 



and 



0.1 < a < 0.2, 



7 > -1 , 



n = 1 + 2a. 



(2.72) 



(2.73) 



In the following sections we will apply these findings to electromagnetic and axionic seeds produced in string cosmology. 
In the axions case we will discuss separately massless and massive perturbations. 



III. SEEDS FROM STRING COSMOLOGY 



In this section we compute the seed functions f p , /„, /„•, and we estimate the Bardeen potentials for electromagnetic 
and axion perturbations, including the case of massive axions. 
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A. Amplification of quantum fluctuations 



We start by recalling the form of the (string- frame) low-energy string effective action |2q| : 

T s eff = J dPxy/\f\e~* (r + g^d^O^ - ^gVg^H^H^ - \g» v g^F^ 



(3.1) 



where we have included the antisymmetric tensor H^a = d[^B va \ and the U{\) gauge field — d^A^. Note that 
this gauge field is typical of what emerges from heterotic string compactification. For gauge fields originating a la 
Kaluza-Klein, the action and the spectra are somewhat different, as discussed in [p7| . 

Upon compactification down to four dimensions, and after introduction of the axion field a by the duality trans- 
formation: 

Hi""* = e~ V^^ct, (3.2) 
one easily arrives at the dimensionally reduced action: 



T s eff = J d^x^Ue-* (r + g^d^d^ - -e^g^d^a ~ \g^ Y° T^) . (3.3) 

The study of tensor (T), scalar-dilaton (SD), electromagnetic (EM) and axion (AX) perturbations is conveniently 
performed defining the external "pump field" , responsible for their amplification. To this aim, we first identify for 
each perturbation the canonical variables ip l , which diagonalize the perturbed action expanded up to second order 
p8|. In a purely metric-dilaton background, such variables are easily found from (3.3) to be: 



= ae-^ 2 h TT = a E h TT , ^ SD = ae^ 2 ^ +..., 

^EM = £ -4>/2 A ^ ^AX = aeH 2 a = aA(J _ (3 4) 

Here h TT denotes the transverse-traceless part of the metric perturbations, the dots in the equation for ^p SD represent 
the additional scalar- metric terms needed to reproduce the gauge-invariant scalar perturbation |28[ , is the scale 
factor in the Einstein frame, and a a in the axion frame [fl2|| . By varying the perturbed action, wc find that the Fourier 
modes ipkiv) °f each of these four perturbations satisfy decoupled, linear equations of the type: 

t>k + U 2 - = , (3.5) 



where Pijf) is the pump field, obtained for each case from eq. (3.4) as: 

p T = pSD^ aE . pE M = e -H2. P AX = aA . (3.6) 
At the beginning of the inflationary era, c hara cterized by an accelerated evolution of the pump field, every per- 



turbation is well inside the horizon and Eq. (3.5) has oscillating solutions, which can be consistently normalized to 
a vacuum fluctuation spectrum. During the whole pre-big bang phase the general solution can be written in terms 
of Hankel functions |29| , with a Bessel index determined by the power that characterizes the background evolution 
(in conformal time) of the pump field. This behaviour has to be matched with the one after the pre-big bang phase 
when, as we assume, the universe becomes radiation-dominated and the dilaton freezes to its present value. In all four 
cases this implies a free Klein-Gordon equation for the canonical variable after the period of accelerated evolution. 
By matching the pre-big bang and radiation solutions of the perturbation equations, we eventually obtain the final 
amplified perturbations during the radiation era. 

For T and SD perturbations the time evolution of the background leads to a spectrum that is in general too 
steep [Q (see also ||o)) to be expected to give any significant contribution to very large scale structures, or to 
temperature anisotropics on the COBE scale. The only way to achieve a reasonable contribution would be to have a 
very long string phase with an almost constant dilaton [^), which is not excluded, in principle, either theoretically or 
phenomenologically, but which looks somewhat unlikely, from both points of view. 

For EM perturbations, however, the situation seems to be more interesting. Consider in fact the transition from 
a pre-big bang phase, with growing dilaton (4> = — 2/3 log |r/|), to the standard radiation-dominated phase with <fi — 
const, and call rji the transition time scale. The electromagnetic fluctuations are directly coupled to the dilaton 
background, in such a way that each polarization mode ipk satisfies at all times, in momentum space and in the 
radiation gauge, the evolution equation: 
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A. = o. 



(3.7) 



In the pre-big bang phase, the general solution of this equation, normalized to a vacuum fluctuation spectrum, can 
be written in terms of Hankel functions of the second kind as: 



^k=»7 1/2 tf< 2) (M)> M = 



V < Vi- 



lli the radiation era we have instead the free plane-wave solution, 

1 



ip k = — [c+(k)e 



V > m ■ 



- 4fc " +c_(fc)e 4fc?? ] 

By matching the two solutions at the transition time rji we easily obtain, for \kr]i\ <C 1 and r\ > 771, 



c± = ±c(fc)e 



Vk 



sin k\rj — rji \ 



\c(k)\ ~ (fc/fci)-^ 1/2 , 



(3.8) 



(3.9) 



(3.10) 



where ki = l/|f?i| represents the maximal amplified frequency (higher-frequency modes are unaffected by the back- 
ground transition) . The associated energy-density distribution of the produced photons is then [pT| : 



dp(k) 
d\og k 



a) M*>l J -lTj (i 



k <ki , fi< 3/2, 



(3.11) 



where /i < 3/2 to avoid photon overproduction which would destroy the homogeneity of the classical background, 
and where the amplitude c(fc) has been estimated modulo numerical factors of order 1. At large times j; > \rji\ we 
thus obtain, in string cosmology, a cosmic background of electromagnetic fluctuations that, for a long enough pre-big 
bang phase with f3 < 2, are characterized by a rather flat spectrum, and could provide the long-sought origin of the 
galactic magnetic fields Jll|| . The amplified fluctuations satisfy stochastic correlation functions, as a consequence of 
their quantum origin. 

Correspondingly, if we consider axionic perturbations, we are led to the canonical equation 



i>k 



k 



OA 
a A 



V>fc = 0, 



(3.12) 



very similar to Eq. (3.7). The same procedure as in the electromagnetic case then leads to the spectrum (3.11) with 



= \r\, where r parametrizes the three-dimensional axion scale factor as a,A(v) ~ rf +1 / 2 . For r = —3/2, in particular 
the axion metric describes a de Sitter inflationary expansion, and the energy density of a massless axion background 
has a flat spectral distribution, dp/d\ogk ~ (ki/a) 4 , as first noted in |l2[ . The value of r depends on the number and 
on the kinematics of the internal dimensions, and the value —3/2 can be obtained, in particular, for a ten-dimensional 
background with special symmetries ]27j . 

In the axion case, however, the low frequency tail of the spectrum is further affected by the radiation — > matter 
transition, as the axion pump field cia is not a constant (unlike the dilaton) in the matte r-dom inated era, where 



a a = a oc 77 . This has important consequences that will be discussed in detail in subsection [II C 



After these preliminary observations we shall now estimate the form of the seed functions for both EM and AX 
seeds. 



B. Electromagnetic seeds 

Here we determine the spectral components of the inhomogeneous stress tensor, for a stochastic background obtained 
by amplifying the quantum EM fluctuations of the vacuum, as discussed in the previous subsection. However, inde- 
pendently of the production mechanism, the results of this section can be applied to any EM fluctuation background 
parametrized by a vector potential that, in momentum space and in the radiation gauge, takes the form 

c(k) 

A,(k,7y) = -^sinfcTy, k t A % = 0, A = . (3.13) 
Ai is a Gaussian random variable which obeys the stochastic average condition: 
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(A i (k)A*(k')) = ^-S 3 (k-k') U 



k 2 



\MKvT- (3-14) 



The above condition has been normalized in such a way that 

£<A j; (k)4*(k')) = {2n)H\k - k>) |A(k,77)| 2 . (3.15) 

i 

Taking into account that the electric component of the stochastic background is rapidly dissipated, because of the 
conductivity of the cosmic plasma J3lj , the seed stress tensor can be expressed in terms of the magnetic field only. 
Setting Bi{k) = ie^ikjA^k), the condition ( |3. 14 ) implies 



(iMk)i?*(k')) = ^V(fc- k') fa - ^) \B(Kv)\\ (3-16) 

where 

|B(k,?7)| 2 = k 2 |A(k,T7)| 2 = /c|c(k)| 2 sin 2 fc?7. (3.17) 

In a process of photon production, the coefficient |c(k)| 2 represents the Bogoliubov coefficient |28| fixing the average 
photon number density, (n(k)}, and is linked to the spectral energy distribution by 

Mkl_ A\ 4 Wfc)) A\ 4 |c(k^ 

dlogfc ~ {a J ir 2 ~ \a) ir 2 ■ (6 
In what follows we shall use for |c(k)| 2 a power-law spectrum, characterized by a cut-off frequency fci, 

|c(k)| 2 = | (Vfcir 2 ^ 1 , , kx, M < 3/2 (3 19) 



This reproduces in particular the spectral distribution (3.11), where /x is fixed by the dilaton growth rate. 

We shall now compute the two-point correlation functions, for the various components of the inhomogeneous stress 
tensor T", associated with the electromagnetic background: 

Q(x, x ') = {t;{x)t;{ x ')) {t;{x)){t;{ x ')) (3.20) 

(no sum over /x, v, and the angular brackets denote stochastic average). The Fourier transform of £ is related to the 
scalar seed functions f p , f v , f v , defined in the previous section. For £q we have, for instance, 

e (x,x>) = J T^e^^ ')|/ p (fc)| 2 . (3.21) 

For Ei = 0, in particular, we have to compute the correlation of a sum of terms that are quadratic in the magnetic 
field. We start considering the energy-density correlation function, 

e (x,x') = (p(x)p(x')) - ({p(x))) 2 , p = -T ° = JH1 (3.22) 

ana 

and compute 

Af^x') = (B 2 (x)B 2 (x>)) - {B 2 ){B 2 ) (3.23) 



where, using the stochastic average ( 3.16 ) and the reality condition B*(k) = B(—k), 



^)) = \l^mk)Hl-%) (3-24) 



2 J (2tt) 3 1 Wl V k 2 

In momentum space, the two-point correlation function for the energy density can be written as a four-point correlation 
function for the stochastic fields (see also fl32"f ). We have, in particular, 
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(B*(x)B 2 (x')) 



d 3 k d 3 k' d 3 p d 3 q i( 

iO V 



k-x+k'-x') 



B i (p)B i (k-p)B j (q)B j (k'-q)). 



(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 
Decomposing the four-point bracket of the Gaussian variables Bj as 
(B l (p)B l (k~p))(B J (q)B ] (k' - q)) + 

+ (B l (p)B 1 ( q ))(B l (k~p)B J (k / ~ q)) + {Bfo)B0 - q)){Bi{k - p)B j (q)), 



(3.25) 



(3.26) 



and using Eq. (3.16), we find that the first term in the above equation is exactly cancelled by the quadratic averages 
(B 2 )(B 2 ), while the other two terms give (no sum over 



A*(x,x') 



d 3 k d 3 p 



* k - Ax 'B(p)| 2 |B(k-p)| 2 U PiPj 



(2tt) 3 (2tt) 3 ' v " P 

where Ax = x — x' . By summing over the vector components we obtain 

A B (x,x')=^A?(x,x') = 



(k-p)j(k-p) 3 

|k- P |2 



(3.27) 



d 3 k d 3 p 



ik-Ax 



iB(p)i 2 |B(k- P )r 



IP ■ (k - P) 
p 2 |k — p| 2 



According to Eq. ( 3.2l| ) , the energy-density spectrum of the electromagnetic seeds is thus determined by 



\f P \ 



d 3 p 

2(87ra 4 ) 2 J (2^ 



|c(p) | 2 |c(k — p)| 2 p|k — p| (l + cos 2 a) sin 2 prj sin 2 |k — p|?y , 



(3.28) 



(3.29) 



where a is the angle between p and k — p. Inserting the power spectrum (3.1£), and defining y = p/fci, z = k/fci, 
the above integral can be written, in polar coordinates, as 



kf 



2(87ra 4 ) 2 (27r) 



dyy 2 ~ 2f * / dxf3~ 2 ^ (l + cos 2 a) sva 2 {ykirj) sin 2 (/3fci7?), 



where we defined x = cosz9, d being the angle between p and k, and 

/3 2 = |z - y| 2 = y 2 + z 1 - 2xyz, cos 2 a = fi~ 2 {y 2 + x 2 z 2 - 2xyz) 



(3.30) 



(3.31) 



The integral of Eq. ( |3.30 ) will be evaluated for \kr)\ = \zktrj\ 1, since we are interested in the large-scale sector of 
the CMB anisotropy, namely in the spectrum of all modes that are still outside the horizon at the time of decoupling. 
For EM seeds these modes give the dominant contribution to the SW effect, as we will see in Section [Fv|. Estimating 
the con tribu tions to the integral from the regions prj <C 1, pr\ ~ 1 and prj 3> 1 , and recalling that fi < 3/2 according 
to Eq. (3.11), we find that the dominant contribution comes from prj S> 1 if ji < 3/4. If 3/4 < fi < 3/2, the integral is 
dominated from its contribution at p ~ k, thus prj < 1 on super- horizon scales. In both cases we obtain for f p a white 
noise spectrum, i.e. |/ P (fc)| 2 ~ constant, but in the second case there is a parametric enhancement (see Appendix B). 
More precisely 



dl{k x laf{k/k x f, 
c 2 (fe 1 / a ) 8 (fc/fc 1 ) 3 (fc 1?7 ) 4 f- 3 , 



fi < 3/4 

3/4 < n < 3/2 , 



(3.32) 



where d p and c p are dimensionless numbers of order 1. Consequently, the energy-density contribution of the EM seeds 
to the Bardeen potentials is, according to Eq. J2.40f) , 



erf 



l/ P |fc 3/2 



4 7 rGd p (a77) 2 (fc 1 /a) 4 (fc/fc 1 ) 3 / 2 , ji < 3/4, 

4 7 rGc p (ar;) 2 (fc 1 /a) 4 (fc/A: 1 ) 3 / 2 (fc 1 ?7) 2 ' 1 - 3 / 2 , 3/4 < yt < 3/2. 



(3.33) 



The contribution of the off-diagonal scalar potential / ff can be similarly obtained by computing the correlation 
function (x, x'), with i ^ j. For purely magnetic seeds, /„ = 0, we find 



fp — 3/p ~ k f-x, 



(3.34) 
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so that the Bardeen potentials, according to Eq. J2.40 ), are always dominated by /„• on super- horizon scales, as 
il 2 f P /U ~ (kr]) 2 < 1- Therefore 

fc 3 / 2 |*-$| ~efc 3 / 2 |M ~ 

47rG^(ar 7 ) 2 (A :i /a) 4 (A : /A :i )- 1 /2(fc l7? )-2 j ^ < 3/4, 

47rGc ;r (ar;) 2 (fc 1 /a) 4 (fc/fci)- 1 / 2 (fci77) 2 ^- 7 / 2 , 3/4 < ^ < 3/2. ' 1 j 

where and are dimensionless numbers of order 1. By assuming that the universe becomes immediately radiation- 
dominated at the physical cut-off scale H\ = ki/ai, such a fluctuation spectrum can be expressed in terms of 
^7 (v) = (Hi / H) 2 (ai / a) 4 , i.e. of the fraction of critical energy density in radiation at a given time 77, and of 
gi = Hi/Mp, the transition scale in units of the Planck mass M p . Denoting with uj = k/a the proper frequency, and 
using p c = 3M 2 H 2 /8ir for the critical density, we obtain 

k ^ (^WW.J-^^/ff)^^, 3/4< M <3/2. (3J6) 



C. Axionic seeds 



As a second example of seed inhomogeneities we will consider a pseudoscalar stochastic background, amplified 
according to the perturbation equation fl3.12| ). 

In the initial, higher-dimensional pre-big bang phase, i.e. for 77 < 771, the solution for the canonical variable ip can 
be written as in Eq. ( |3.8| ), with fi = \r\ < 3/2, as discussed previously. In the radiation era, i.e. for 771 < 77 < r] eq , 
the effective potential cia/cia is vanishing, as <p =const and a <~ 77, and ip is given by the plane-wave solution (3.9). In 
the final matter-dominated era, i.e. for 77 > rj eq , we have a ~ rj 2 , and 0,4/0,4 = 2/?7 2 . The plane- wave solution is still 
valid for modes with k > k eq — r)~ q , that are unaffected by the last transitions. Modes with k < k eq feel instead the 
effect of the potential in the matter era, and the general solution of Eq. ( |3.12| ), for those modes, can be written as 



w ,) = |(m|;)im}i) 



l [(A + B)J 3/2 - i(A - B)Y 3 / 2 ] , k<k eq , r,>r hq . (3.37) 



Vk 

Here J 3 / 2 and Y3/2 are Bessel functions of argument kr] (we follow the convention s of P9| ). 

The matching of the solutions at r\\ determines the coefficients c±(k) as in eq. ( 3.10| )7The matching at rj eq gives 

J 4 + B~c(k)(fe ?e ,) -1 , A - B ~ c(k) (k Veq ) 2 , (3.38) 

so that the contribution of J 3 / 2 to ip^ is always dominant with respect to the Y3/2 contribution, both for kr) > 1 
and kr) < 1. In the matter-dominated era, i.e. for 77 > ?7 eg , we can thus approximate the produced stochastic axion 
background as follows: 

c(k) 

er(k, 77) ~ — — sin fc?7, k > k eq , 



i*Jk 

c(k) / fc 

c(k) / k 

a\fk V feg 



-1 

\2 



(kr)) , fc < k^, kr] < 1, 

k < k eg , kr] > 1. (3.39) 



The correlation functions for the various components of the stress tensor, 

T; = d^a - ~8» (d a a) 2 (3.40) 
can be computed by exploiting the stochastic average conditions of the Gaussian variables a, & and <jj = djcr, 
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<a(k)a*(k')) = (2vr)^ 3 (fc - fc'JEiCM), 

(<j(k)<7*(k')) = (27r) 3 5 3 (fc-fc')S2(k,r7), 

(tr i (k)o-;(k')) = k i k j (2Tr) 3 5 s (k-k , )X 1 (k,r l ), 

(^■(k^k')) = -(<r(k)<T*(k')) = ifypTr) 3 ^* - fc')S 3 (k,77) 



where, according to Eq. ( |3.39 ) 



£i(M) 



|c(k)| 2 ^ 



lc(k)l 2 f k ^ ' 



ka 2 \ k, 



k < k eq , krj > 1 



(3.41) 



(3.42) 



£ 2 (k,77)~fc 
~0, 



|c(k) 



|c(k)| 2 / /,• * 



k > k e q , 

k <c k^q^ kt] *\ 1) 
< fc eg , krj > 1, 



(3.43) 



£ 3 (M) 



k(k)r 

a 2 

0, 

k(k)r 



k > k, 



eq 



eq i 

< fc e g, /cry < 1, 
^ kgq . kfj ^ > 1 . 



(3.44) 



Following the same procedure as the one used for EM seeds, and collecting all contributions to the correlation function 
of the axion energy density, 



1 



Pa = 2^ ^ + ' 
we obtain from £q(x, x') that the energy density spectrum is determined by 



(3.45) 



k 3 \M 2 ( - 



2fc 3 



d 3 p 



(2a 2 ) 2 J (2tt) 3 
2p.(k-p)S 3 (p)S3(k-p) 



S 2 (p)S 2 (k - p) + |p • (k - p)| 2 S 1 (p)E 1 (k - p) 



(3.46) 



In order to evaluate this integral outside the horizon, in the region krj < 1, we must distinguish two cases, [i < 3/4 
and fi > 3/4. In both cases, by separate integration in the ranges < p < rj^ 1 , rf x < p < k eq , k eq < p < k±, we find 



a white noise spectrum, \f p \ ~ const. In particular (see Appendix B): 

*.3/2 |f | - J ^(fci/a) 4 (fcAi) 3/2 [l + ^(Wfci) 2 (^) 2 ' I+1/2 ] 

lp 'Uy " l ^(fci/a) 4 (fcAi) 3/2 (fce g /fci) 2 (fcir/) 2 ^ +1 / 2 , 



H < 3/4 
3/4 < /j < 3/2, 



(3.47) 



where c p ,d p ,5 p are dimensionless numbers of order 1. The same power spectrum is also obtained for the scalar 
velocity potential f v , associated to the axion seeds. An explicit computation gives in fact kf v ~ krjf p so that the 
contribution of f p and f v to the Bardeen potential are both of the same order, namely: 



af 



\f p \k^~er,^\f v \k^ = 



_ J 4 7 rG^(ary) 2 (fc 1 / a ) 4 (fc/fci) 3/2 [l + ^(* efl /*i) 2 (A 1 »?) 2 '* +1/a ] , 
4 7 rGc-M 2 (fc 1 /a) 4 (fc/A :i ) 3 / 2 (fc eg /fc 1 ) 2 (fc 1 r;)^+ 1 /2, 



fi < 3/4 

3/4 < /j < 3/2. 



(3.48) 
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We will now consider the anisotropic stress potential f w , defined according to (2.23) by: 



v 4 /.= 



2M^' 



CTiCTj 



-Sijidka) 



Sin Bid j 



(3.49) 



Summing all contributions in the two point correlation function, we find 



d 3 k 



2M 4 J (2tt) 3 



d 3 p 
(2^)3 



ik- Ax 7,4 



1 



1 



p |k — p| Ei(p)Ei(k — p) ( cos iJcos 7 — - cos cos 7 cos a + - cos a), 



(3.50) 



where a and 7 are, respectively, the angles between p and k, p and k — p an d k and k — p. The integral over p is 
of the same type as t he integral for the energy density spectrum (see Eq. ( ft.46] )), and gives for k 2 f^ the same white 
noise spectrum (3.47) as for f p (modulo numbers of order 1) , since 



k m \Uk)\{^j M 2 -k^\f p (k)\{^ 



(3.51) 



On super-horizon scales the contribution of / w to the Bardeen potentials is always dominant with respect to the f p 
contribution since, from the above equation, 



V 2 f P ~ (kvTU- 



(3.52) 



In the whole range kr\ < 1 we can thus estimate the scalar perturbation spectrum, induced by massless axion seeds, 
through the contribution to the Bardeen potentials. We find 



4 7 rG<(ary) 2 (fc 1 /a) 4 (fc/A :i )- 1 /2(fc 1?? )-2 [1 + ^(^/^a^^+i/a] > ^ < 3/4j 



4 7 rG<(a7 7 ) 2 (fc 1 /a) 4 (fc/fci)- 1/2 (fce 9 /fci) 2 (fci7 7 ) 2 ^ 3 / 2 , 



3/4 < n < 3/2, 



(3.53) 



gln^iuo/^y^i^/H)-* [1 + ^M) W# ) 2Al+1/2 ] , I* < 3/4, . . 

5 2 fi 7 (ry)(^ 1 )- 1 /2( Weg / Wl )2 (LJl/i/) 2p-3/2 j 3/4</x<3/2. 



where d^, 5° are dimensionless numbers of order 1. As we will see in Section [V, the dominant contribution to the 
SW effect now comes, for each mode, from the time of reentry 7/ ~ 1/k. 
Let us finally discuss the case of massive axions, with 

T; = d^d v u - ~5; [{d a af - m 2 a 2 } , (3.55) 



and a primordial distribution again characterized by the index /i, as in Eq. ( |3.19 ). The mass term directly contribute 



to f p and f p , and only indirectly to the off-diagonal potentials f v , f^. We are interested in the axion perturbations 
that may be relevant to the large-scale CMB anisotropy, namely in the modes that are outside the horizon at the 
decoupling era, k < aHd ec - If, for these modes, the mass contribution is negligible, ma < k < aHdec, then the AX 
seed functions and the corresponding Bardeen potentials are the same as in the massless case (see before). We will 
thus concentrate our discussion on the case in which the axion mass is large enough, so that all modes outside the 
horizon at the equilibrium epoch are already non-relativistic: 

ma > aH eq > k. (3.56) 

In this case we may neglect the effects of an additional axion production in the matter-dominated era, since a 2 m 2 > a/a 
at r\ > rj eq . The axion fluctuations are amplified by the inflation — * radiation transition, but are to be evaluated in 
the non-relativistic regime (77 > r] eq ), where the mass contribution is already important. 

For non-relativistic, super-horizon modes, the Fourier components of the axion field become (see the non-trivial 
calculation reported in Appendix C): 

c(k) / k \ 1 ^ 2 / Hi \ 4 ^ 4 . fm\ f m ^ ' 
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where the initial distribution c(k) is still given by Eq. ( 3.19| ). Here k m is the limiting frequency re-entering the horizon 
at the same time as it becomes non-relativistic, i.e. k m /a m — H m — m. Indeed, we are assuming that at the transition 
scale Hi the mass term is completely negligible, m <C Hi, and all modes are relativistic. As the proper momentum 
is red-shifted, the mo des b ecome non-relativistic whe n m = k/a = to, and re-enter the horizon when H = u>. 

For the axion field (4.14) the stochastic conditions (3.41) arc still valid, but the squared amplitude (3.42), averaged 
over time scales m/H 3> 1, now become 



Sl( P ,7?) = 



|c(k)l* 
2ma 3 



Hi 

m 



1/2 



-2-^S 2 (p,?;) = — £ 3 (p,7?). 

i a ma 



(3.58) 



For the case we are considering, the contribution of / w to the Bardeen potentials is always negligible with respect 



to rf f p . An explicit computation gives, in fact, 



rffp/fn ~ m/H > 1, 



(3.59) 



where the last inequality is a consequence of ( |4 . 1 3| ) . In addition, the mass contribution to the AX energy density 
dominates with respect to the momentum contribution, since m > k/a. The energy density correlation function thus 
becomes: 



(as |<r(fc)| = ma\a(k)\), and gives, using Eq. (4.15) 



$(x,x') = m 4 ((a 2 (x)a 2 (x')) - (a 2 (x)) 2 ) 
d 3 p 



(3.60) 



Exfo^k-p) 



mHi ( ki 
8tt 2 (~a~ 



dyy 



2-2/1 



dx(3 



-2/i. 



(3.61) 



[II B 



where x, y and (3 are defined in Section 

It should be noted that the above expression for the spectrum is only valid if /i > 3/4. Only in that case, in fact, 
is the integral over ^dominated by the contribution of the lower boundary, p/ki — > 0, and is the use of Eq. ( 3.58j ) for 
the axion spectrum appropriate. In the opposite case, we have to take into account the different spectrum of non- 
relativistic sub- horizon modes, for p > k m , and possibly of relativistic modes in the high-frequency limit p — > ki. In 
both cases we obtain, for /i < 3/4, a white noise spectrum and a negligible contribution to the large-scale anisotropy, 
as we will see in the next section. 

We will thus concentrate on the case 3/4 < fi < 3/2. For this case the integral (3.61) is estimated in Appendix B, 
and we obtain 



k'lUn-) =<mHi[^) (- 



a I 



ki 



6-4/i 



3/4 < (i < 3/2 , 



(3.62) 



where c m is a dimensionless number of order 1. The corresponding Bardeen spectrum is: 



erf 



\f P \k 3/2 



AicGv^UnirUnlh) 1 ' 2 [ — ) [ 



3-2/i 



9l c m 



n«r(w). 



(3.63) 



We may note that ^H 2 evolves in time like a 3 , so that, during the matter-dominated era (when H 2 oc a 3 ), the 
Bardeen potential \& remains frozen at the value reached at the time r\ eq of matter-radiation equilibrium. Using 
{H 1 /H eq ) 2 (a 1 /a eq ) 3 = (a eq /ai) = (H 1 /H eq ) 1 / 2 , wc obtain for i] > 7j eq , 



|<l>| . f-^)' J f 



3-2/1 



(3.64) 



The CMB anisotropy induced by the EM and AX seeds discussed here will be analysed in the next section. 
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IV. CMB FLUCTUATIONS FROM PRE-BIG BANG SEEDS 



For electromagnetic seeds, with the assumption that the electric field is already dissipated away at recombination, 
we find that the seeds are generically suppressed by a factor (hr/ dec) 2 , and the anisotropic stress f v dominates over the 
density contribution f p (see the discussion at the end of Section II A ). By contrast, for massless axionic perturbations, 
there is no (h]d ec ) 2 suppression for f p , while there is one for f n . For large wave numbers which enter the horizon before 
matter and radiation equality, EM and AX seeds lead to similar amplitudes. Consequently, if the convolution leading 
to fn is dominated by s mall scale contributions, fi < 3/4, the two cases give similar geometric scalar perturbations 
through Eq. ( |2.40| ). 

However, on large scales, krj eq < 1, the additional axion production during the matter-dominated era leads to an 
enhancement by the factor {T)/T) eq ) ■ This changes the time-dependence of the Bardeen potentials and has important 
consequences as we will see below. 



A. Electromagnetic seeds 



The scalar metric perturbation spectrum induced by EM seeds is reproduced in Eqs. ( [3.35] ) and (3.36). Comparing 
with our parametrization in terms of a and 7 (see Eqs. ( |2.59| ), ( 2.60| )) we find 



-4, M < 3/4 

2/i-ll/2, 3/4<^<3/2 { > 

7/2, [i < 3/4 

5-2/x, 3/4<^<3/2 K ' 



and 



^= c -(£) ( fc iW 2 (4-3) 

in both cases /i < 3/4, /i > 3/4 (modulo numbers of order 1). 

Since 7 + 1 < 0, in both cases the seeds decay fast enough outside the horizon, and our analysis of Section [n] applies. 
However, in both cases 7 + a — —0.5, which leads to the spectral index n = 2 , i.e. to a spectrum that grows too fast 
with frequency to fit the results of COBE observations, see Eqs. ( |2.71 ), ( 2.72 ). 



The quadrupole amplitude is given by Qrms-ps = \J (5/47r)C2To, which has been measured |33| to be Qrms-ps — 
(18± 2)/xAT. This leads to 

C 2 = (1.09 ± 0.23) x lO" 10 . (4.4) 



From Eq. (2.7C), using a + 7 = — 1/2, kiij eq = (Hi/ 'H eq ) x l 2 , gi = Hi/M p , and setting I = 2, we obtain: 



62 ~ 10(47r) 4 (7 + 1) 2 [hTJ {-^) {-^J ■ (45) 

Compatibility with the COBE normalization, C2 ^ 10 -10 , thus implies 

(6 - a) log 10 g x < 55(a - 2) - 6 + log 10 ( 7 + l) 2 - log 10 c 2 (4.6) 

(we have used H eq /M p ~ 10 -55 , and rjdec ~ Veq ~ 10 _2 ?7o)- This important constraint is easily satisfied by a growing 
seed spectrum, fx < 3/2, i.e. a > 2. In the limiting (and most unfavorable) case fj, = 3/2, a — 2, 7 = —5/2, the 
above condition reduces to 

log 10 9l < - 1.4 - 0.5 log 10 (V . (4.7) 

Even in this limiting case there are no stringent constraints on the typical inflation scale of the "minimal" pre-big 
bang scenario @J^,§, expected to approach the string mass scale M s as g\ = H\/M p ~ M s /M p . Indeed, the limiting 



condition (4.7) is marginally compatible even with the maximal expected value Hi ~ M s , since [j34| 

10" 2 < M s /M p < lQ" 1 . (4.8) 

To conclude, the EM fluctuations seem to lead to a scalar perturbation spectrum that grows too fast with frequency 
to contribute in a significant way to the observed large-scale anisotropy. The positive aspect of our result is that there 
are no significant constraints from the COBE normalization to the production of seeds for galactic magnetic fields, 
which remains allowed as discussed in MM. 
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B. Axionic seeds 



Let us first consider massless axions. If fj, < 3/4, the situation is like in the electromagnetic case. The CMB fluctu- 
ations induced have the wrong spectrum, but their amplitude is sufficiently low to avoid conflict with observations. 

If3/4</t<3/2 the situation becomes radically different. Comparing Eq. (3.53) with the ansatz ( 2.59] ), ( p. 60 ) we 
obtain, due to the additional factor (r]/r] eq ) 2 , 



7 = 2/i - 7/2, 

a = -7- 1/2 = -2/i 



3 . 



(4.9) 
(4.10) 



In the limiting case \i = 3/2 this yiel ds 7 = —1/2 and a — 0, which corresponds to a Harrison-Zel'dovich spectrum of 
CMB fluctuations, according to Eq. (2.73), with an amplitude 



9l 



(4.11) 



(we have absorbed into g\ all dimensionless numerical coefficients of order one appearing in the spectrum ( 3.53 )). 
Note that f p leads to a Bardeen potential with the same a, but with 7 = 2/i — 3/2. However, since again 7 > — 1, the 
contribution to the SW effect is the same for f p and / w (see Section ||). 

The normalization of the axion spectrum to the COBE amplitude (4.4), according to Eq. ( 2.65| ), imposes the 
condition 



G 



SW 



' M 2 (Mo) 



-2a 



4 / w 

91 u 



6-4^ 



10 



-in 



which implies 



logio 5i 



164- 116/i 

1 + 2/Lt 



(4.12) 



(4.13) 



(again we have absorbed numerical coefficients into gi, and we have used ujq ~ 10 -18 Hz 
to 0,0|). On the other hand, the allowed range for the spectral index (see Eq. (2.73), 
fx < 3/2 (required to prevent over-critical axion production), implies 

1.4 < //< 1.5. 



uji ~ g{ 10 11 Hz, according 
combined with the condition 



The combination of ( |4.13| ), ( fL14| ) leads to 

3 x 10~ 3 < g x = {Hx/Mp) < 2.6, 



(4.14) 



(4.15) 



which is perfectly compatible with the identification Hi ~ M s (see Eq, (|4 



A stochastic background of massless axions, produced in the context of the pre-big bang scenario, is thus a possible 
viable candidate for a consistent explanation of the large-scale anisotropy observed by COBE. The important difference 
between AX and EM seeds is the non-conformal coupling of the axions to the metric, that leads to an additional 
amplification of perturbations after the matter-radiation equality. 

Another interesting case is that of a massive axion background, for which the /„- contribution to the Bardeen 
potentials is negligible when the super- horizon modes are already non-relativistic at the time of decoupling, m > Hd ec - 
As seen in the previous section, one then obtains constant Bardeen potentials, with 7 = 0, a = (n — l)/2 = 3 — 1\x 
and 



<A 



1/2 



(4.16) 



(see Eq. ( |3.64 ), where we have set c m = 1). A flat Harrison-Zel'dovich spectrum is again possible in the limiting case 
fx = 3/2. The amplitude of perturbations, however, is enhanced by the factor (m/ ' H^ ) 1 / 2 , so that the value of the 
axion mass has to be bounded, to avoid conflicting with the COBE normalization (4.4). 

The allowed range for the spectral index, 0.8 < n < 1.4, an d the condition /i < 3/2 to avoid over-critical axion 
production, again imply for the parameter /i the allowed range (4.14). In addition, the present axion energy density 
is constrained by the critical density bound, f2 CT (?7o) < 1, imposed at the peak frequency u) m of non-relativistic modes 
(see Appendix C) . Actually, an even stronger condition is required for the validity of our perturbative approach, which 
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neglects the back-reaction of the axionic seeds on the expansion of the universe. Using Eq. (C25) we thus impose the 
condition 



ri CT (w m ,77o) 



9\ 



H,, n 



1/2 



m 
Hi 



2-/1 



< 0.1, 



which implies 



5 55 

(2 - (i) [log 10 {m/H eq ) - log 10 gi - 55] + - log 10 g x + — < -1. 



(4.17) 



(4.18) 



In order to find a possible AX mass window compatible with the COBE data, we now impose the normalization 
Ci ~ 10~ 10 on the massive axion spectrum Eq. ( |3.64 ). From Eq. ( p. 65 ) we obtain 



H 



eq 



6-4^ 



10" 



from which 



y. ~ [164 - log 10 (m/fl- e ,) - 41og 10 9l ] /116 



(4.19) 



(4.20) 



(we have used uj\/ojq ~ 10 29 , neglecting the weak dependence of u>i on the tran sition scale gi). By eliminating 
fi in terms of m and gi, according to the above equatio n, th e constraints ( 4.14 ) and ( 4.18 ), plus the condition 
to > Hdec ~ H eq (assumed for the validity of the spectrum (3.64)), determine an allowed region in the plane (to, Hi) 
as follows: 



10- 10 {M p /Hxf < m/H eq < lO 16 {M v /H x f , 



to > H„ 



(4.21) 



[68 + log 10 (m/7J eg ) + 41og 10 5i] [log 10 (m/H eq ) - 55- log 10 £fi] + 58 (55 + 5 log 10 g x ) < - 1. 
For a typical string cosmology scale, Hi ~ M s ~ (10 _1 — 10 _2 )M p , we thus obtain the maximal allowed window: 



10" 27 eV <, m <, 10~ 17 eV. 



(4.22) 



As illustrated in Fig. 1, the window is shifted towards higher values of mass as the final inflation scale is lowered, and 
the seed condition (4.17) becomes important only at low inflation scales, H\/M p < 10 -7 . The stringent upper limit 
we obtained for the mass can be traced back to the simplest model of b ackg ro und u sed in this paper, that gives the 
same slope for the axion spectrum at low and high frequency (see eqs. (024), ( |C26D ). It is not excluded that higher 
values of the mass may become possible in a more complicated model of background, giving a steeper high frequency 
spectrum. 



V. CONCLUSIONS 



In this paper we have considered the possibility that, in a string cosmology context, the large-scale temperature 
anisotropics may arise from the contribution of seeds to the metric fluctuations, and not from the direct amplification 
of the metric fluctuations themselves, as in the conventional inflationary scenario. We have discussed, in particular, 
two cases: one in which the seeds are EM vacuum fluctuations amplified by the growth of the dilaton field, and one 
in which the seeds are AX vacuum fluctuations amplified by the time evolution of a higher-dimensional background. 

In the case of EM perturbations we have found that the induced angular power spectrum of AT/T grows too fast 
to be compatible with COBE observations. However, the contribution of the seeds to the large-scale anisotropy may 
be consistently imposed to be negligible, without constraining in a significant way the basic parameters of the pre-big 
bang models. 

Massless AX perturbations, unlike EM perturbations, are also affected by the radiation — > matter transitions. This 
changes the time dependence of the seed contribution to the Bardeen potentials and, due to the integrated Sachs- 
Wolfe effect, a flat or slightly tilted blue spectrum of temperature anisotropies can be induced, compatible with present 
COBE observations. Scale-invariant massless axion seeds thus appear as possible promising candidates for structure 
formation. Determining in more details the CMB anisotropy spectrum also on smaller angular scales requires however 
numerical simulations, which we defer to a future research project. 

For massive AX seeds the situation is qualitatively different if the mass is such that all modes outside the horizon at 
the time of decoupling are already non-relativistic. In that case the contribution to AT/T is controlled by the axion 
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log io ( m/H eq) 

FIG. 1. The phenomenologically allowed region is to the left of the curve fi CT = 0.1, to the right of the vertical dashed line 
m — H eq , and lies within the full lines n = 1, n = 1.4, to avoid conflicting with present COBE observations (n < 1 is excluded 
by over-critical axion production). The shaded area defines the allowed mass window for an inflation scale Hi = M a , typical of 
string cosmology. 



mass, and a slightly tilted blue spectrum is still compatible with the amplitude and the slope measured by COBE, 
provided the axion mass is inside an appropriate window, in the ultra-light mass region. Higher values of masses may 
become possible in models with more complicated backgrounds. 

At smaller angular scales, an axionic origin of CMB anisotropics should lead to acoustic peaks in the spectrum, 
with a structure different from that of the standard inflationary scenario. This may in principle allow a test of models 
with axionic seeds through the very accurate observations of the CMB anisotropy planned in the near future jj5|. A 
thorough investigation of this possibility is postponed to a future paper. 
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APPENDIX A: SACHS- WOLFE COEFFICIENTS FOR POWER-LAW SPECTRA 



Assume that the Bardeen potentials are given by power-law spectra as in Eq. (2J 

* - * = { c^f' * > i ' c[k) = ^ k ~ 3/2{ k,klT ' 

where x = krj, xq — fcryo, Xdec — krj dec . The SW contribution to the angular coefficients Ci is given by 



C, 



SW 



= A . ; ,_ i kl dk,_ 



2n 



where 



= je{x - 1) + / x~<j' e (x - x)dx , 



(Al) 



(A2) 



(A3) 



and a prime stands for the derivative of je with respect to its argument. 

We concentrate here on the case where 7 + 1 < 0. Furthermore, we are interested in the situation where the integral 
in Eq. (A2) is dominated by large scales (small values of fc), and therefore Xd ec <C 1. In that case the integral I{k) is 
dominated by its value at the lower bound: 
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7+1 



21+ I 



■ ( \ 1+1 ■ I \ 

Ji-i{xa) - 2j—3i+i{xa) 



This leads to the following expression for the CVs 

2(7+1) 

>W _ I 'Idee • 

J l 

71 \ 



2 / %ec\ 1 /, w2a / da; 2 

~ TT, — ia&iVa) / ^o 



rf^O 2(q+7+1) 



(2£+l) 2 



if-i(^o) 



2f (£ + 1) . 



|l + 7l 2 * 

1)2 



2(7+1) 



(2^ + 1)5 

(Mo) -1 



(l+l) 2 
(2£+l) 



X 



p ad _ m+A rm + (^ +1 ) 2 jo) 

(2£+l) 2 * (2£ + l) 2 £ (2^+l) 2 ' 



where, setting jj> = we find (Ref. pq|, number 6.574) for a + 7 < 0, 



r(U _ 



7T 

2 JO 
7T 



dxa; 2(Q+7) j|_ 1/2 (a;) 

r(-2(a + 7 ))ro? + a + 7 ) 



1 D ~T 



2 2- 2 («+^)[r(-(a + 7) + l/2)] 2 r(£ - (a + 7)) ' 

poo 

:; ; ^ 2(a+ ^_ 1/2 (^+3 /2 (z) 

z Jo 

7T 



r(-2(a + 7 ))r(£+l + a + 7 ) 



r(3) _ 



2 2- 2 ( Q +7)r(-(a + 7) - l/2)r(3/2 - (a + j))T(l + 1 - (a + 7)) ' 

/>oo 

— / dxa; 2| ' Q+7 ' 1 j| +3 ^ 2 (a;) 







7T 



r(-2(a + 7))r(^ + 2 + a + 7) 



2 2- 2 («+7)[r(-(a + 7) + l/2)] 2 r(£ + 2 - (a + 7)) ' 
Finally, combining the above results, we obtain the result given in Eq. ( [2.70 ) 



pjSW 



N 2 



r(-2(a + 7 )) 



Vdec 

m 

I 2 e-a-j 2i(£+l) 1/2 + a + 7 



2-2(a+ 7 )( 7 + 1)2 y(1/2 - a - 7 ) 2 



2(7+1) 



-2a r(l+l + a + 7 ) 

r^ + la-7) 



'+1) 2 ^+l + a + 7 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



_(2£+ l) 2 ^ + a + 7 (2£+l) 2 l/2-a- 7 (21 + l) 2 £ + 1 - a - 7 

f 2) 

It is interesting to note that, for a + 7 = —1/2, the mixed term vanishes, which is indeed what happens 
case of electromagnetic seeds (see Section |Tv| ). 



(A9) 
in the 



APPENDIX B: THE SEED FUNCTIONS 



1. Electromagnetic Seeds 

For purely magnetic seeds, all t he se ed functions can be approximately determined by the energy density correlation 
function £g, which leads to Eq. (3.3C). The contribution of super-horizon modes (kr) < C 1) to the spectrum can be 
estimated in the limit z = k/k\ — > 0. In this limit (3 — > y, cos 2 a — > 1, and the integral (3.30) reduces to 



dy y 2 -^ sin 4 (y k lV ), ft < 3/2. 



(Bl) 



The dominant region of integration is easily shown to be y ~ 1 for /1 < 3/4 and ykirj ~ 1 for 3/4 < /1 < 3/2. This 
gives 

(fci/a) 8 (fc/fci) 3 , M < 3/4 

(fc 1 /a) 8 (A;/A;i) 3 (ibi?7) 4 ' t - a , 3/4 < fi < 3/2 ' 



(B2) 



modulo numerical factors of order one. This coincides with the result reported in Eq. (3.32) 
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2. Massless Axions 



For massless axions, the seed spectral functions are determined by the integrals (3.46), (3.50). The various terms 



appearing in the integrands turn out to give comparable contributions, so let us concentrate on the typical term 

7=^ |d 3 pp 2 |k-p| 2 I] 1 (p)S 1 (k-p). (B3) 

We distinguish different integration regions: 0<p<k, k<p< r\~ x , r\~ x < p < k eq , k eq < p < k\. The dominant 
integration regions depend on the value of \i but, for all /i < 3/2, they always lie at p > rj^ 1 > k. This is the reason 
why we always obtain a white noise spectrum. On the other hand, the behaviour in rj depends on which region of p 
dominates. Specifically we find: 

1) F or 3/4 < < 3/2 the leading contribution to I comes from p ~ 77" 1 , and gives the single term appearing in eq. 



1} For /i < 3/4 the leading contribution comes either from p ~ k\ (giving the first term in the square brackets of 
&), or (for li very close to 3/4) from p — 77 1 (giving the second term in the same brackets). 



3. Massive Axions 



For massive actions, the energy density spectrum is determined by Eq. ( 3.61 ) , with 3/4 < /i < 3/2. This integral is 
dominated by the region p ~ k, and its rough behaviour can be easily obtained this way. For a more precise evaluation 
we proceed as follows: the angular integration gives 

^ (t)' - Teroi (£)' (£)' jf * k. - »r» - c + *r»i ■ M 

Defining t = y/z we obtain 



where, after some manipulation |36| 



POO 

A= dt t 1 ' 2 " [(1 - t) 2 - 2 " - (1 + tf 
Jo 



o4/j — 4 

-^^r(2 - 2/i)r(2/i - 3/4) [cos2tt(/^ - 1) - 1] (B6) 



and 



fl= / dt t 1 - 2 ^ [(1 - i) 2 " 2 ' 1 - (1 + tf- 2fl ] . (B7) 

Jl/z 

By evaluating this second integral in the limit z —> 0, we obtain 

B ~ z 4 ^ 3 < A. (B8) 

so that 

= 16*^1) U) UJ ^ 3/4<,<3/2, (B9) 



as reported in eq. (3.62). Note that there is no singularity for \i = 1, as 



r(2 2//) 2 ^ _ _ = _ 1)2 = const (B10) 

At-i (m-i) (a* -i) 2 
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APPENDIX C: NON-RELATIVISTIC CORRECTIONS TO THE AXION SPECTRUM 



For a massive-axion perturbation a, the string frame action 

3 = \ J d A x^—ge* [{d,af - mV 2 ] , 

in a conformally flat background, can be written in terms of the canonical variable 

tp = za, 



<t>/2 



as 



S = - I d 6 xdr) 



ae 



i? - {d^f + V - m 2 a 2 ^ 2 

z 



(CI) 



(C2) 



(C3) 



(the dot denotes differentiation with respect to the conformal time 77) . The Fourier modes ipk satisfy the perturbation 
equation 



k h mV } ip k = 0. 



(C4) 



We shall consider the background transition at rj = r\\ from an initial pre-big bang phase in which the axion is 
massless, to a final radiation-dominated phase in which the dilat on f reezes to its present value, and the axion acquires 
a small (in string units) mass. For r\ > r\\ the solution of Eq. (C4) depends on the kinematics of the pump field z 
and, after normalization to an initial vacuum spectrum, it can be written in terms of the second-kindHankel functions 



Mv)=V 1/2 H (2) (kr ] ). (C5) 
In the radiation era, r\ > rji, the "effective potential" z/z is vanishing, and the perturbation equation reduces to 

4> k + (k 2 + aV) </>fc = 0, (C6) 

where we have put 

a = mH 1 al, (C7) 



„,2 2 2 2 

ma = q ■ 1) 1 , 



using the time behaviour of the scale factor, a ^ rj. 

Assuming that the mass term is negligible at the transition scale, m <C k/a, we can match the solution (C5) to the 
plane-wave solution 



[c+(k)e 



-ikr] 



c_(fc)e 



ikr/l 



and obtain: 



c± 



±c(k)e 



\c(k)\ ~ {k/h 



n-M-1/2 



(C8) 



(C9) 



(We are neglecting, for simplicity, numerical factors of order 1, which are not very significant in view of the many 
approximations performed. Their contribution will be included into an overall numcricalcoefficient in front of the final 
spectrum.) In the rclativistic regime, the amplified axion perturbation then takes the form: 



v(Kv) = -^7= sm (fc?7), 



a\fk 



(CIO) 



used in Section III C for the massless- axion case. 



In the radiation era the proper momentum is red-shifted with respect to the rest mass, and all a xion modes tend to 
become non-relativistic. When the mass term is no longer negligible, the general solution of Eq. (C6) can be written 
in terms of parabolic cylinder functions ]2l| . For an approximate estimate of the axion field in the non-relativistic 
regime, however, it is convenient to distinguish two cases, depending on whether a mode k becomes non-relativistic 
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inside or outside the horizon. Defining as k m the limiting comoving frequency of a mode that becomes non-relativistic 
(km = fna rn ) at the time it re-enters the horizon (k m = H m a m ) 7 we find, in the radiation era, 



We will thus consider the two cases k ^> k m and k <C k n 



m \ 

#i7 



1/2 



(Cll) 



In the first case, we rewrite the perturbation equation (C6) as 



d 2 ^ k 
dx 2 



— -b}i/) k = Q, 



x = rj 



(2«) 1/2 , 



b = k 2 /2a, 



and we give the general solution in the form 

ip = AW(b,x) + BW(b,-x) 



(C12) 



(C13) 



where W(b, x) are the Weber parabolic cylin der fu nctio ns (see | p9| , chap. 19). In order to fix the integration constants 
A and B we shall match the solutions ( C13 ) and ( |C10| ) in the relativistic limit 
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m 2 a 2 a 2 r] 2 



» 1. 



(C14) 



In this limit, as we are considering modes that become non-relativistic when they are already inside the horizon, 



__) ~L„(_ 6 )»i, 

a 



(C15) 



we can expand the W functions for b large with x moderate j29|. Matching to the plane- wave solution ( CIO ), we 
obtain A = 0, and 



In the opposite, non-relativistic limit x 2 3> \4b\, the expansion of the Weber functions gives |2£ 



c(k) 



(ar)) 



■<) . (m\ 



(we have used x 2 /4 = marj/2 ~ m/H). The corresponding axion field is (inside the horizon) 



o-(k,?7) 



c 



(k) . 

, sm 



ita (j?) 



(C16) 



(C17) 



(C18) 



Consider now the case of a mode that becomes non-relativistic when it is still outside the horizon, k <C /c m . I n thi s 
case, we cannot use the large \b\ expansion as \b\ < 1, and it is convenient to express the general solution of Eq. ( C12 ) 

as 



ip = Ayi(b,x) + By 2 (b,x) , 

where y\ and yi are the even and odd parts of the parabolic cylinder functions 
relativistic limit x — ► 0, gives A = and 

/ k \ 1/2 

■0fc~c(k)f— j y 2 (b,x). 

In the non-relativistic limit £ 2 ^> |6| we use the relation 0] 

zj 2 ~ [W(6, z) - W(6, -a:)] ~ 4= sin X' 

4 



(C19) 

Matching to jCK| ), in the 
(C20) 

(C21) 



which leads to 
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to 



c(k) 



2\ V4 



sin — 



(C22) 



Using Eqs. flC15| ) and (^Tl]) for /c 2 /a, we finally arrive at the non-relativistic axion field presented in Eq. (|4.14|): 



cr(k, 7]) 



k 



c(k) ( 



1/2 



•Hi 



1/4 



sm 



n 



(C23) 



For later use, it is also convenient to define the spectral energy density in critical units, Q a (u>) = d(p/p c )/dbxu), 
associated with the stochastic axion backgr ound in the three different regimes defined before. 
For relativistic modes we find, from Eq. (CIO), 



3-2^ 



Hi 
H 



/ai 
V a 



m < uj < u>x. 



For modes that becomes non-relativistic after re-entry we find, from Eq. ( |C18[) , 



Q. a {uj) ~ g\ 



Hi \LOl 



Hi 
H 



UJm < UJ < 771. 



For modes that becomes non-relativistic before re-entry we find, from Eq. (C23), 



fi ff (u>) ~ g\ 



rn 
H~i 



1/2 / x 3-2u 
' UJ x 



Hi 
i? 



UJ < UJ ri 



(C24) 



(C25) 



(C26) 



The last two spectral distributions are constant durin g the matter-dominated era, and the last one corresponds to the 
spectrum of the Bardeen potentials, as given in Eq. (3.63). 



[1] G. Veneziano, Phys. Lett. B 265, 287 (1991); M. Gasperini and G. Veneziano, Astropart. Phys. 1, 317 (1993); Mod. Phys. 
Lett. A 8, 3701 (1993); Phys. Rev. D 50, 2519 ( 1994). An updated collection of papers on the pre-big bang scenario is 



available at 



http : //www. to . inf n. it/~gasperin/ 



[2] R. Brandenberger and C. Vafa, Nucl. Phys. B 316, 391 (1989); G. Veneziano, Ref. 0; A. A. Tseytlin, Mod. Phys. Lett. 
A 6, 1721 (1991); K. A. Meissner and G. Veneziano, Phys. Lett. B 267, 33 (1991); Mod. Phys. Lett. A 6, 3397 (1991); A. 
Sen, Phys. Lett. B 271, 295 (1991) ; S. F. Hassan and A. Sen, Nucl. Phys. B 375, 103 (1992); A. A. Tseytlin and C. Vafa, 
Nucl. Phys. B 372, 443 (1992); M. Gasperini and G. Veneziano, Phys. Lett. B 277, 256 (1992). 
[3] G. Veneziano, Phys. Lett. B 406, 297 (1997); A. Buonanno et al., Phys. Rev. D 57, 2543 (1998). 

[4] R. Brustein and G. Veneziano, Phys. Lett. B 329, 429 (1994); N. Kaloper, R. Madden and K. A. Olive, Nucl. Phys. B 
452, 677 (1995); Phys. Lett. B 371, 34 (1996); R. Easther, K. Maeda and D. Wands, Phys. Rev. D 53, 4247 (1996); R. 
Easther and K. Maeda, Phys. Rev. D 54, 7252 (1996); M. Gasperini, J. Maharana and G. Veneziano, Nucl. Phys. B 472, 
349 (1996); S. J. Rey, Phys. Rev. Lett. 77, 1929 (1996); M. Gasperini and G. Veneziano, Phys. Lett. B 387, 715 (1996); 
M. Gasperini, M. Maggiore and G. Veneziano, Nucl. Phys. B 494, 315 (1997); R. Brustein and R. Madden, Phys. Lett. B 
410, 110 (1997); Phys. Rev. D 57, 712 (1988). 

[5] M. S. Turner and E. J. Weinberg, Phys. Rev D 56, 4604 (1 997); M. Maggiore and R. Sturani, Phys. Lett. B 415, 335 
(1997); N. Kaloper, A. Linde and R. Bousso, |hep-th/9801073| . 

[6] G. F. Smoot et al, Ap. J. 396, LI (1992); C. L. Bennet et al., Ap. J. 430, 423 (1994). 

[7] M. Gasperini and M. Giovannini, Phys. Lett. B 282, 36 (1991); Phys. Rev. D 47, 1519 (1993); M. Gasperini and G. 
Veneziano, Ref. fl; R. Brustein, M. Gasperini, M. Giovannini, V. F. Mukhanov and G. Veneziano, Phys. Rev. D 51, 6744 
(1995). 

[8] R. Brustein, M. Gasperini, M. Giovannini and G. Veneziano, Phys. Lett. B 361, 45 (1995); R. Brustein, M. Gasperini and 

G. Veneziano, Phys. Rev. D 55, 3882 (1997). 
[9] M. Gasperini, Phys. Rev. D 56, 4815 (1997). 
[10] R. Durrer, Phys. Rev. D 42, 2533 (1990). 

[11] M. Gasperini, M. Giovannini and G. Veneziano, Phys. Rev. Lett. 75, 3796 (1995); D. Lemoine and M. Lemoine, Phys. 
Rev. D 52, 1955 (1995). 



26 



[12] E. J. Copeland, R. Easther and D. Wands, Phys. Rev. D 56, 874 (1997); E. J. Copeland, J. E. Lidsey and D. Wands, 

Nucl. Phys. B 56, 506 (1997). 
[13] M. Gasperini, M. Giovannini and G. Veneziano, Phys. Rev. D 52, 6651 (1995). 
[14] R. Durrer and M. Sakellariadou, Phys. Rev. D 56, 4480 (1997). 
[15] R. K. Sachs and A. M. Wolfe, Ap. J. 147, 73 (1967). 
[16] R. Durrer, Fund, of Cosmic Physics 15, 209 (1994). 
[17] J. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[18] H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1(1984). 

[19] R. Durrer and N. Straumann, Helv. Phys. Acta 61, 1027 (1988). 

[20] J. M. Steward and M. Walker, Proc. R. Soc. London A 341, 49 (1974). 

[21] R. Durrer, A. Gangui and M. Sakellariadou, Phys. Rev. Lett. 76, 579 (1996). 

[22] R. Bond and G. Efstathiou, MNRAS 226, 655 (1987). 

[23] P. J. E. Peebles, The Large Scale Structure of the Universe (Princeton University Press, Princeton, 1980). 

[24] G. F. Smoot and D. Scott, in L. Montanet et al., Phys. Rev. D 50, 1173 (1994) (update 1996). 

[25] H. R. Harrison, Phys. Rev. D 1, 2726 (1970); Y. B. Zel'dovich, Mon. Not. Roy. Ast. Soc. 160, 1 (1972). 

[26] R. R. Metsaev and A. A. Tseytlin, Nucl. Phys. B 293, 385 (1987). 

[27] A. Buonanno ct al., JHEP01, 004 (1998). 

[28] V. F. Mukhanov, A. Feldman and R. Brandenberger, Phys. Rep. 215, 203 (1992). 

[29] M. Abramowitz and I. A. Stegun, Handbook of mathematical functions (Dover, New York, 1972). 

[30] M. Gasperini, Relic dilatons in string cosmology, in "Proc. of the 12th Int. Conf. on General Relativity and Gravitational 

Physics", Rome, 1996, editors M. Bassan et al. (World Scientific, Singapore, 1997), p. 181. 
[31] M. S. Turner and L. M. Widrow, Phys. Rev. D 37, 2743 (1988). 

[32] M. Giovannini, Teoria quantistica delle perturbazioni in cosmologia di stringa, PhD Thesis, Turin University (January 
1996). 

[33] A. J. Banday et al., Astrophys. J. 475, 393 (1997). 
[34] V. Kaplunowsky, Phys. Rev. Lett. 55, 1036 (1985). 

[35] M. Bersanelli et al., COBRAS/SAMBA, Report on the phase A study, ESA Document D/SCI(96) 3 (1996). 
[36] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (Academic Press, New York, 1965). 



27 



